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ABSTRACT 


c 


This  tedinical  report  develops  a  theoretical  bit  error  rate  expression  for  a  Kaiman  filtering 
technique  that  is  used  for  filtering  narrowband  interferers  out  of  direct  sequence  spread 
spectrum  signals.  The  approach  is  based  on  the  digital  phase-locked  loop  Kalman  filter 
and  is  close  to  optimum  in  so  far  as  demodulating  an  FM-type  of  interferer.  Because  the 
interference  is  assumed  to  be  much  stronger  than  either  the  signal  or  noise,  the  Kalman 
filter  locks  onto  the  interference  and  produces  estimates  of  its  phase  and  envelope.  Several 
assumptions  are  made  in  the  development  of  the  theoretical  bit  error  rate  expression 
and  these  are  expanded  upon.  It  is  shown  and  illustrated  that  theory  agrees  well  with 
simulation  as  long  as  the  bandwidth  of  the  low  pass  filter  used  to  determine  the  envelope 
of  the  interference  is  much  less  than  the  sampling  rate,  otherwise  the  interference  estimate 
at  the  input  to  the  interference  canceller  will  have  significant  cross- correlation  with  the 
signal,  leading  to  some  cancellation  of  the  spread  spe'-trum  signal. 

R^SUMfe 


Ce  rapport  elabore  une  expression  mathematique  donnant  le  taux  d’erreur  de  bits  quand 
un  filtrage  de  Kalman  est  utilisee  pour  la  suppression  d’interferences  sur  des  signaux 
a  spectre  etale  par  sequence  directe.  L’approche  repose  sur  un  filtre  de  Kalman  as- 
servi  numeriquement  par  verrouillage  de  phase  et  s’avere  quasi-optimum  quant  a  la 
demodulation  d’une  interference  de  typ>e  MF.  Puisqu’on  pr^uppose  que  I’interference 
est  plus  forte  le  signal  ou  que  le  bruit,  le  filtre  de  Kalman  se  verrouille  sur  I’interference 
et  permet  I’estimation  de  sa  phase  et  de  son  enveloppe.  On  elabore  sur  la  nature  des 
hypotheses  qui  doivent  etre  posees  lors  du  developpement  de  I’expression  gouvernant  le 
taux  d’erreur  de  bits.  On  demontre  eloquemment  que  la  thwrie  est  en  parfaut  accord 
avec  les  simulations  tant  que  la  largeur  de  bande  du  filtre  passe-bas,  servant  a  circon- 
scrire  I’enveloppe  de  I’interference,  est  inferieure  au  taux  d’echantillonnage.  Autrement, 
I’estimation  de  I’interference,  a  I’entrw  meme  ou  elle  doit  etre  cancellw,  sera  affect^ 
de  fa^on  importante  ptir  des  intercorrelations  avec  le  signal,  amenant  un  evanouissement 


partiel  du  signal  a  spectre  etale. 


EXECUTIVE  SUMMARY 


This  technical  report  develops  a  theoretical  bit  error  rate  expression  for  a  Kalman  fil¬ 
tering  technique  that  is  used  for  filtering  narrowband  interferers  out  of  direct  sequence 
spread  spectrum  signals.  These  signals  are  used  extensively  in  military  communication 
systems.  The  technique  described  herein  applies  equally  to  both  Electronic  Support  Mea¬ 
sures  (ESM)  systems  and  direct  sequence  spread  spectrum  communication  syterns.  In  the 
former  application,  the  ESM  system  may  be  attempting  to  intercept  the  spread  spectrum 
signal,  but  the  narrowband  interference  may  be  hampering  this  effort.  In  the  latter  appli¬ 
cation,  the  spread  spectrum  communication  system  may  require  additional  assistanr**  to 
suppress  ».he  interference.  Since  the  open  literature  has  been  devoted  to  this  latter  case, 
the  material  presented  here  focuses  on  the  communications  application. 

One  of  the  attributes  of  direct  sequence  spread  spectrum  communication  systems 
is  their  ability  to  combat  interference  or  intentional  jamming  by  virtue  of  the  system's 
processing  gain  inherent  in  the  spreading  and  despreading  process.  The  interference  can 
be  attenuated  by  a  factor  as  high  as  this  processing  gain.  In  some  cases,  however,  the  gain 
is  insufficient  to  effectively  suppress  the  interferer,  leading  to  a  significant  degr£idation  in 
communications  manifested  by  an  increase  in  bit  error  rate.  If  the  ratio  of  interference 
bandwidth  to  spread  spectrum  bandwidth  is  small,  the  interference  can  be  filtered  out 
to  enhance  system  performance.  However,  this  is  at  the  expense  of  introducing  some 
distortion  onto  the  signal.  This  process  of  filtering  is  sometimes  referred  to  as  interference 
excision. 

The  K2dman  filtering  approach  to  excision  is  based  on  the  digital  phaise-locked 
loop  Kalman  filter  which  is  close  to  optimum  in  a  minimum  mean-squared  error  sense  for 
demodulating  an  FM  interferer.  Because  the  interference  is  assumed  to  be  much  stronger 
than  either  the  signal  or  noise,  the  Kalman  filter  locks  onto  the  interference  and  produces 
an  estimate  of  the  phase  and  envelope  of  the  interference. 

Several  assumptions  are  made  in  the  development  of  the  theoretical  bit  error 
rate  expression  and  these  are  expanded  upon.  It  is  shown  mathematically  and  illustrated 
that  theory  agrees  well  with  simulation  as  long  as  the  bandwidth  of  the  low  pass  filter 
used  to  determine  the  envelope  of  the  interference  is  much  less  than  the  sampling  rate.  If 
this  condition  is  not  met,  then  the  interference  estimate  at  the  input  to  the  interference 
canceller  will  have  significant  cross-correlation  with  the  signal,  leading  to  a  certain  degree 
of  signal  cancellation. 
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1.0  INTRODUCTION 


The  problem  of  interference  suppression  in  direct  sequence  spread  spectrum  communica¬ 
tions  sy terns  is  well-known  (2,  3,  4,  5]. 

A  suppression  technique  based  on  the  Kalman  filter  was  presented  in  [6j.  There, 
a  function  related  to  the  interference  is  assumed  to  have  been  generated  by  a  second 
order  nonlinear  state-space  system  to  which  the  extended  Kalman  filter  equations  are 
applied.  This  results  in  what  has  been  termed  the  digital  phase-locked  loop  (DPLL) 
(1,  7).  This  approach  to  interference  suppression  is  close  to  optimum  for  the  suppression 
of  constant  envelope  or  FM-types  of  interferers  and,  therefore,  can  potentially  provide 
better  performance  than  the  recursive  least  squares  (RLS)  algorithm  which  was  analyzed 
for  the  swept  tone  in  [5,  8]. 

The  objective  of  this  technical  report  is  to  derive  a  theoretical  equation  describing 
the  bit  error  rate  performance  of  the  DPLL  and  to  compare  this  theoretical  expression 
to  simulated  results.  The  interference  is  assumed  to  be  of  the  FM-type  generated  by  the 
state-space  model  defined  in  [6,  1]. 

The  outline  of  this  technical  report  is  as  follows.  Section  2.0  describes  the  spread 
spectrum  communications  model.  Section  3.0  details  the  interference  model  used  in  gen¬ 
erating  the  interference  and  the  Kalman  algorithm.  Section  4.0  presents  the  interference 
estimator.  Section  5.0  develops  the  theoretical  bit  error  rate  equation.  This  expression 
is  compared  to  simulated  results.  Finally,  Section  6.0  concludes  the  technical  report, 
suggesting  areas  for  further  resejurch. 

2.0  COMMUNICATIONS  MODEL 

The  basic  elements  of  the  BPSK  PN  spread  spectrum  receiving  system  are  shown  in 
Fig.  1.  The  received  waveform  r{t),  consisting  of  a  spread  spectrum  signal,  additive 
white  Gaussian  noise,  and  narrowband  interference  is  applied  to  a  bandpass  filter  with 
the  transfer  function  Hhp{f),  whose  output  is  defined  as 

u(t)  =  s(t) -b  n(t) -1- t(t).  (1) 

The  bandpass  filter  Hbp{f),  for  the  application  considered  here,  is  assumed  to  be  a  filter 
matched  to  a  chip  and  centered  at  the  angular  carrier  frequency  a?o  of  the  spread  spectrum 
signal.  The  spread  spectrum  signal  is  defined  as 


s(t)  =  a{t)  cos(u>ot) 


(2) 


rn 

Correlator 


Figure  1:  Spread  spectrum  communications  model. 

where 

=  (3) 

k 

In  Elq.  (3),  D*  is  a  sequence  of  data  bits  of  amplitude  (±1)  and  duration  7),  seconds  (bit 
rate  jR*  =  l/T*),  and  bk{i  —  kTt,)  is  the  PN  sequence  pattern  for  Ihe  bit,  i.e., 

L 

f>k{t)  =  ^Ckjq{t-jTc)  (4) 

with  L  being  the  number  of  pseudo  random  chips  per  bit,  or  the  processing  gain,  Ckj  is 
the  code  sequence  for  the  bit,  q{t)  representing  the  basic  chip  pulse  of  energy  Ec,  and  Tc 
the  reciprocal  of  the  chip  rate  Rc. 

The  noise  n(<)  is  Gaussian  and  has  a  power  spectral  density 

Sn(/)  =  w) 

where  No/2  is  the  power  spectral  density  of  the  assumed  white  Gaussian  noise  from  the 
channel.  The  band  of  interference  is  defined  as 

i{t)  =  I(t)cos(iLfoi  +  0(0)  (6) 

where  1(0  is  the  interference  amplitude  and  0(0  is  the  phase  modulation.  It  has  been 
assumed  that  the  effect  of  the  bandpass  filter  Hbp(f)  is  negligible  on  the  interference  *(t). 
This  would  not  be  true  in  reality,  since  the  bandpass  filter  will  induce  some  amplitude 
modulation  on  the  constant  envelope  interference.  The  rate  of  change  of  the  envelope 
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would  therefore  depend  on  the  rate  of  change  of  the  instantjuaeous  frequency,  aiid  the 
amplitude  deviation  would  depend  on  the  frequency  offset  from  wo  and  the  interference 
bandwidth.  The  fact  that  this  report  considers  the  amplitude  to  be  of  constant  envelope 
a-  ‘he  output  of  the  bandpass  filter  corresponds  to  the  best  case  situation. 

Referring  to  Fig.  1,  the  output  u(t)  of  the  bandpass  filter  is  bandpass 

sampled  and  applied  to  a  limiter/bandpass  filter  and  interference  ^timator. 

Consider  the  bandpass  sampler  first.  The  analog  signal  u(f)  from  Eq.  (1)  is 
sampled  at  /,  =  2Rc  (mf,  =  u»o/25r  +  Rcf2  for  some  integer  m).  The  resultant  sampled 
signal  is,  therefore, 

Ufi  =  -Sfi  -b  n„  -f-  in,  (7) 

where  s„  consists  of  the  sequence  {. . .  ,0,(-l)’*a„,0,  where  the  are 

of  energy  Ee  and  coded  according  to  CkjDk  for  the  chip  in  the  transmitted  bit. 
Tin  ^  are  uncorrelated  Gaussian  noise  samples  of  variance  =  Ec{No/2),  and  is  the 
sampled  version  of  Eq,  (6).  The  samples  u„  are  applied  to  the  interference  estimator  and 
interference  canceller. 

The  interference  estimator  produces  an  estimate  t„  of  the  interference  which  is 
removed  from  the  sampled  input  u„.  The  output  of  the  summer  is  decimated  and  sign- 
reversed,  resulting  in  a  baseband  error  signal  e„.  This  error  signal  is  correlated  with  the 
PN  sequence.  The  output  of  the  correlator  is  integrated  and  bit-detected. 


3.0  INTERFERENCE  MODEL  AND  KALMAN  ALGORITHM 


This  section  presents  the  state-space  models  of  the  interference.  The  development  here 
follows  that  in  [l], 

A  block  diagram  of  the  analog  state  space  model  which  generates  the  analog 
version  of  a  function  related  to  in  Fig.  1  is  shown  in  Fig.  2;  this  function  is  represented 
by  bp{ty.  The  state-space  equations  are  defined  as 


Xi{t) 

X2{t) 

= 

’  0 

0 

d ' 

—a 

Xi{t) 

X2(<) 

+ 

0 

y/KZ 

u;{t). 


(8) 


where  a:i(f)  and  xj(f)  are  the  phase  function  and  modulating  signal  state  variables  of  the 


^Coherent  bandpass  sampling  has  been  assumed,  so  that  the  in-phase  component  is 
(ni,n/v/5)coe(nx/2)  and  the  quadrature  component  is  (nj  n/v^)  8in(nT/2) 
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»<t) 


Figure  2:  Block  diagram  of  the  state  space  model. 


interference,  respectively.  In  compact  form  Eq.  (8)  becomes 

x{t)  =  Ax{i) gw{t).  (9) 

The  interferer’s  modulating  signal  X2{t)  is  assumed  to  be  generated  by  applying  Gaussian 
white  noise,  variance  to  an  amplifier  with  gain  y/KZ  followed  by  a  first  order 

low  pass  filter  of  bandwidth  a  =  2icaj.  The  steady  state  variance  of  the  coloured  noise 
at  the  output  of  the  lowpass  filter  is  unity  if  K^,  =  2a.  The  output  is  multiplied  by  the 
frequency  deviation  constant,  d,  and  then  integrated  to  yield  the  phase  function 

Xi(t)  =  d  f  X2iT)dT.  (10) 

JQ 

An  FM  type  of  signal  (interferer)  is  generated  by  phase  modulating  a  carrier,  u/q,  which 
has  been  set  equal  to  the  carrier  of  the  spread  spectrum  signal,  since  any  offsets  from  u-'o 
can  be  incorporated  into  Xi(0i  yielding 

y(t)  =  v/2cos(u;o<  +  (11) 

This  signal  is  corrupted  by  white  Gaussian  noise  v{t)  with  power  spectral  density  Noi,al2, 
generally  different  from  No/2  defined  earlier  for  the  communications  model  in  Fig.  1.  The 
result  is  am  FM  signal  defined  as 

bp{ty  =  y{t)  +  v{t).  (12) 

Given  the  observation  process  bp'{t)  from  Eq.  (12),  the  objective  is  to  estimate 
the  state  xi(t)  and  the  related  state  X2(t),  as  well  as  y{t)  in  Eq.  (11).  One  approximate 
technique  reserved  for  nonlinear  estimation  problems  [9]  is  the  extended  Kalman  filter 
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algorithm. 

Before  the  algorithm  is  presented,  the  discrete  forms  of  Eqs.  (8)  and  (9)  will  be 
stated,  namely, 

*n+-l  =  d>*n  +  (13) 

where  n  represents  discrete  time, 


’  1  ^(1  -  e-“^*)  ' 

f  '  " 

1 _ 

0 

[0  J 

(14) 


and 


to. 


=  / 

JnT. 


(n^-l)r. 


<t>gv}{z)  dz. 


(15) 


In  Ekj.  (14),  ^  =  dja  and  is  termed  the  bandwidth  expansion  ratio  in  units  of  volts~^  [10] 
and  Tf  is  the  sampling  interval.  In  Eq.  (15),  to„  is  a  stationary  zero-mean  white  Gaussian 
vector  sequence  whose  covariance,  (i.e.,  £{to„io^}  where  the  superscript  t  refers  to  the 
matrix  transpose),  is  given  by 


Vw  -  <p99*<t>du. 

JnT, 


Finally,  Eq.  (12)  becomes, 


+  On, 

where  y„  is  the  sampled  version  of  Eq.  (11),  i.e.. 


(16) 


(17) 


y„  =  \/2cos(wonr,  -f  Xi,„). 


(18) 


It  has  been  assumed  that  bj/{t)  from  Eq.  (12)  has  been  applied  to  a  bandpass  filter  of 
bandwidth  B,  whose  output  has  been  sampled  at  a  rate  so  as  to  yield  uncorrelated  noise 
samples  Vn  of  variance  al  =  [7].  It  should  be  noted  that  6p[,  in  Eq.  (17)  is 

closely  related  to  5pn  in  Fig.  1.  The  difference  is  due  to  the  fact  that  in  5p„,  as  shown 
later,  the  noise  is  not  additive  as  in  Eq.  (17),  but  is  contained  in  the  argument  of  the 
cosine  function  of  y„  in  Eq.  (18)  due  to  the  hardlimiting  process.  For  large  interferen  v 
to-signal-plus-noise  conditions,  however,  the  effect  of  the  additive  noise  in  Eq.  (17)  will 
be  similar;  thus  this  model  is  valid. 

The  noise  covariance  matrix,  V^,  is  obtained  by  substituting  Eq.  (14)  and  g 
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from  Eq.  (9)  into  Eq.  (16)  and  integrating  each  term,  yielding 


^  [iirh  -  3  + 

2a  (l  -  1  - 


(19) 


where  7  =  2w/{aTt),  which  is  the  ratio  of  the  sampling  rate  /,  =  1/T,  to  the  interferer’s 
modulating  signal  bandwidth  a/.  The  discrete  form  of  the  interference  model  is  shown  in 
Fig.  3. 


Figure  3:  Discrete  form  of  the  interference  model. 


The  demodulator  in  Fig.  3  is  the  Kalmem  filter  algorithm  listed  in  Table  1  [6,  1]. 
In  the  algorithm  in  is  defined  as  the  estimate  of  the  true  state  vector  in  is  defined  as 
the  post  error  Xn  —  *nj  is  the  Kalman  gain  vector,  Vf,  =  E{inin]  is  the  post  error 
covariance  and  V is  the  prior  error  covariance.  Furthermore, 
/njn-i  1®  Itnown  as  the  “coupled  equation”,  since  the  estimate  of  the  interferer’s  phase  is 
contained  in  it.  However,  this  equation  can  be  further  approximated  as  noted  in  [1].  The 
approximation  is 


/-I 


I 


V- 


(20) 


The  block  diagram  of  this  demodulator  is  illustrated  in  Fig.  4  (also  referred  to  as 
the  digital  phase-locked  loop  or  DPLL).  It  produces  an  estimate  of  the  phase  = 

^n|n-i  ^d  a  signal  in  quadrature  to  y„  in  Eq.  (18);  this  signal,  defined  as 

V  COn  =  -y/2  am{uJonTt  +  K\n-i ),  (21 ) 


is  analogous  to  the  voltage-controlled  oscillator  (VCO)  signal  of  a  phase-locked  loop. 

Performance  curves  of  the  demodulator  can  be  obtained  fl]  using  the  baseband 
simulation  model  with  unit  sampling  interval,  i.e.,  T,  =  1  or  /,  =  1  Hz.  A  few  examples 
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Table  1;  The  Kalman  Filter  Algorithm. 


Initial  Conditions: 


*_i  -  E{x.i}  = 


mi 

mj 


=£{*.iiL,} 


Do  for  n  =  0  to  T: 

^nln-l  =  [  1  0  ]  <ftXn-l 

^_i  _ _ 1  -  cos(2fa;on  +  2gn|n>i) _ 

“  cos(2wbn  +  2^„|n-i) 


Vi.  =  Vi...  .  -  Vi 


*»|n— 1 


1  0 
0  0 


Kn  =  -Vi. 


l/<rl 

0 


V2sin(u;bn  +  0n|n-i) 


*„  =  <l>Xn-i  +  Kr,  -  V2coa{uJon  + 


Figure  4:  State  space  representation  of  the  DPLL. 
will  now  be  presented. 

The  theoretical  performance  curves  were  obtained  using  the  prior  error  covari¬ 
ance  equation  in  Table  1  (i.e.,  The  experimental  curves  were  obtained  using  the 

model  in  Fig.  3  for  a  range  of  values  for  d  and  a/  was  fixed  at  0.0003125  Hz.  The 
results  are  illustrated  in  Fi^.  5  and  6  and  refer  to  the  filter  form  of  the  Kalman  algorithm 
which  is  listed  in  Table  1.  The  modified  Kalman  algorithm,  which  is  similar  (except  for 
the  definitions  in  ^  and  V,„),  is  described  in  [6]. 

The  results  are  illustrated  in  Figs.  5  and  6  and  indicate  reasonable  agreement 
between  theory  and  simulation.  It  must  be  pointed  out  that  the  theory  does  not  predict 
the  threshold  at  which  the  variance  of  the  phase-noise  increases  dramatically;  this  must 
be  determined  from  experimentation.  As  cited  in  [1],  an  approximate  threshold  is  a 
phase  variance  of  0.25  radians^ /second*,  (i.e.,  1/Vi„  „  =  4).  The  results  for  the  modified 
Kalman  algorithm  as  shown  in  Fig.  6  exhibit  a  higher  threshold.  Otherwise,  from  the 
perspective  of  phase-noise  variance,  the  two  algorithms  perform  almost  identically. 

4.0  INTERFERENCE  ESTIMATOR 

Referring  to  Fig.  1,  consider  the  branch  containing  the  limiter.  Here,  Un  is  applied  to  a 
limiter/bandpass  filter.  The  input  to  the  limiter  referenced  to  the  interference  is  redefined 

as  _ 

«n  =  y/[In  +  -I-  [n',,,,  +  o'j.J*  cos(won  +  -|-  (22) 


8 


9 


JNR  (=  2/Ta/^o4.) 


Figure  6;  Performance  curves  using  the  modified  Kalman  filter  in  the  filter  mode. 


0 


where 


^it,n  =  arctaxi 


/  ”3.n/v^  +  Q2.n  \ 


(23) 


is  a  Qoise-Uke  phase  fluctuation  on  the  interferer’s  phase  and  is  due  to  the  noise  and 
spread  spectrum  signal.  The  terras  n[  ^,  a[  ^,  and  a,  „  are  in-phase  and  quadrature 
components  of  the  noise  and  spread  spectrum  signal  with  respect  to  the  interference  phase 
9n,  i.e., 


”i,n  =  «i.n  coe(^«)  +  nj-n  sin((9n)  (24) 

-  ^3.n  C06(5„)  -  n,,„  8in(dn)  (25) 

<*i.«  =  anCos(^«)  (26) 

a;,,  =  -n„8in(d„).  (27) 


If  the  interference  is  of  constant  envelope  and  /  »  (sn  -f-  n„),  then  Eq.  (23)  simplifies  to 


^i*,n  ^ 


The  output  of  the  limiter/bandpass  filter  is  [11] 


4A' 

6p„  = - cos(u;on  + 


(28) 


(29) 


where  A'  is  the  limiter’s  output  level.  This  signal  is  redefined  as  (letting  A*  =  \/2x/4  for 
convenience) 

bpn  =  V^coslwon  +  (30) 

It  should  be  noted  that  for  large  interference-to-noise  ratios  in  which  the  interference  is 
of  constant  envelope,  ^u,n  in  E<1-  (30)  is  approximately  Gaussian  [11].  The  sampled  signal 
in  Eq.  (30)  is  what  is  processed  by  the  Kalman  filter,  which  estimates  the  phase  Bn  of  the 
interference.  The  phase  estimate  is  denoted  as  ^nin-i  • 

The  interference  estimator  shown  in  Fig.  1  is  detailed  in  Fig.  7.  The  Kalman 
algorithm  from  Fig.  4  produces  the  VCO  sign2il  defined  in  Eq.  (21)  and  shifted  by  x/2 
radians.  This  signal  is 

yn  =  V2  cos(<non  -|-  9n\n-i )  (31) 

and  can  be  used  as  a  basis  for  estimating  a  sampled  version  of  the  envelope  of  i{t)  (i.e., 
In  in  Eq.  (6)).  The  output  of  the  first  multiplier  is  a  baseband  term  and  is,  using  Eq.  (7), 
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Figure  7:  Block  diagram  of  the  interference  estimator, 
excluding  the  \/2  factor, 

K  =  /nCOs(d„  - 

+  Hn  C08(won  +  ^n|n-l )  +  <*«  COs(5„|„_i ) 

+  o„cos(2ci;ori +  (32) 

Equation  (32)  consists  of  four  terms:  the  first  term  is  related  to  the  desired  envelope  of 

the  interference;  the  second  is  approximately  Gaussian  baseband  noise  [10];  and  the  third 
and  fourth  terms  are  noise-like  terms  emanating  from  the  spread  spectrum  signal.  The 
fourth  term,  because  of  the  sampling  rate  conditions  discussed  in  [1],  is  essentially  filtered 
out  by  the  low  pass  filter  of  bandwidth  Blpf  <  0.50  Hz  and,  therefore,  will  be  ignored  in 
the  baseband  simulations  to  be  discussed  in  the  next  section.  The  term  /^/V^  is  filtered, 
resulting  in  the  estimate  of  the  interference  envelope,  /„/ y/2.  Combining  this  with  in 
Eq.  (31)  and  shown  in  Fig.  7  yields  the  estimate  of  the  interference 

*  A  A 

*f»  ~  InCOs(l^n  -f-  dn|n— l))  (30) 

which  is  subtract^  from  Un  as  illustrated  in  Fig.  1. 

5.0  THEORETICAL  BIT  ERROR  RATE  PERFORMANCE 

For  the  development  of  the  theoretical  bit  error  rate  expression,  one  can  refer  to  Figs.  1 
and  7.  In  Fig.  1,  the  error  signal  e„  prior  to  despreading  (but  following  decimation  and 
sign  reversal)  will  be  used.  This  error  is  defined  as 

A 

A 

*  3n 
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+  Hn  +  Al„. 


(34) 


In  Eq.  (34)  the  signal  Sn  has  been  set  equal  to  a„  which  was  defined  following  Eq.  (7),  and 
Hn  has  been  set  equal  to  the  in-phase-noise  tenn  ni,n<  It  is  assumed  that  the  interferer  is 
of  constant  envelope  and  that  the  baseband  representation  of  it  is  used,  i.e., 

i„  =  /co8(0„).  (35) 


Furthermore,  it  is  generated  using  the  interference  model  in  Fig  3. 

From  Fig.  7,  the  interference  estimate  »„  is 

in  =  I  COB0n\n-l)  (36) 

where  _ 

/  =  ^/cOs(£„)  +  \/Ecl2nB,n  +  Sn  COs((?„|„_i  )j  *  htPF.n-  (37) 

InEq  (37), 

Cn  —  ~  ^n|n-l  (38) 

ftB.n  is  a  baseband  noise  term  defined  as 

nfl,n  =  ni,„  cos(^„|„_i )  +  n2,„  8in(^„|„_i )  (39) 

where  the  terms  ni,n  and  nj.n  are  the  in-phase  and  quadrature  noise  terms  of  prior  to 
decimation.  Furthermore,  hipp.n  is  the  impulse  response  of  the  low  pass  filter  in  Fig.  7 
and  ♦  is  the  convolution  operator.  With  Ec  of  unit  energy  and  Eq.  (37)  substituted  into 
Eq.  (36),  the  interference  estimate  becomes 

tn  =  {  [/  C08(€„)]  ♦  htppn}  cos(0„j„_i )  + 

+  {  [nfl,„/\/2  +  S„  COs(5„|„_i)]  ♦  htPF,n}  cos(^„|n-i).  (40) 

This  interference  estimate  consists  of  two  terms:  (a)  the  first  term  is  closely  related  to 
the  interference  in  Ek}.  (35);  (b)  the  second  term  is  low  pass  filtered  noise  and  signal, 
modulated  by  an  estimate  of  the  phase  variations  of  the  interference. 

The  first  term,  because  it  is  a  function  of  ,  which  in  turn  is  a  function  of 
the  past  data,  is  uncorrelated  with  the  signal  and  noise,  s„  and  Un-  The  same  cannot  be 
said  about  the  second  term,  since  it  contains  noise  and  signal  terms  at  time  n.  However, 
the  degree  of  cross-correlation  will  be  a  function  of  the  bandwidth  of  the  interference 


13 


characteristics  and  the  low  pass  filter  with  impulse  response  hipp^n- 

To  illustrate  the  point,  consider  a  simple  example  which  is  close  to  the  worst 
case  situation  in  which  the  low  pass  filter  in  Fig.  7  is  of  a  rectangular  shape  of  bandwidth 
Blpf  =  0.50  Hz.  Furthermore,  assume  the  interference  has  a  very  naurow  bandwidth, 
close  to  a  single  tone,  offset  from  the  spread  spectrum  signal  by  a  small  amount  Sf 
Hz.  If,  for  illustrative  purposes,  the  phase-noise  is  considered  to  be  negligible,  then 
*  2x6 fn,  so  that  the  signal  term  in  Eq.  (32),  anCO8  0„|„_i,  after  low  pass  filtering, 
will  have  a  truncated  spectrum  offset  from  0  Hz  by  ±6 /  Hz.  This  new  signal,  when  mixed 
by  y„  in  Fig.  7,  will  result  in  a  shifting  back  to  baseband  of  the  truncated  spectrum,  with 
significant  signal  energy  in  the  neighbourhood  of  0  Hz,  i.e.,  in  the  region 

-Blpf  -{-Sf  <  f  <  Blpf  -  Sf- 


After  decimation  this  baseband  term  becomes 


A  A 

{[fln  COs(d„(„_i)]  ♦  UlPf])  COs(5„|n_i) 


where,  as  noted  earlier,  =  s„.  One  can  see  that  the  amount  of  signal  energy  in  the 
region  around  0  Hz  diminishes  as  Blpf  decreases  from  its  maximum  value  of  0.50  Hz  and 
as  6/  increases.  A  similar  situation  occurs  with  the  noise. 

Therefore,  consider  the  case  where  the  bandwidth  Blpf  is  sufficiently  small  so 
that  the  second  term  in  Eq.  (40)  is  weakly  correlated  with  Sn  and  n„.  With  these  simpli¬ 
fications,  the  signal-to-noise  ratio  at  the  output  of  the  despreaider  can  be  approximated 
by 


SNRo 


LNo/2  +  LE{Ail} 
L 

No/2  +  } 


(41) 


where  the  derivation  of  Eq.  (41)  is  similar  to  the  derivation  of  Eq.  (59)  in  [5]  for  the  linear 
prediction  case.  In  deriving  Eq.  (41),  it  has  been  assumed  that  Ai’n  is  weakly  correlated 
with  s„  and  Un  for  the  reasons  stated  above,  that  s„  is  uncorrelated  with  n„,  and  that  the 
chips  of  the  PN  code  over  a  bit  interval  are  uncorrelated.  Finally,  because  the  energy  per 
chip  was  set  to  unity,  and  the  number  of  chips  per  bit  is  L,  the  energy  per  bit  is  Ei,  =  L. 
Thus,  for  no  residual  interference  Ai,  Eq.  (41)  reduces  to  the  standard  equation  for  BPSK 
signals  in  AWGN. 

If  it  is  assumed  that  the  noise  plus  residual  interference  have  Gaussian  statistics 
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[12],  the  bit  error  rate  can  be  represented  by 

Pi  =  ^trfcyfsNRji,  (42) 

where  erfc{x)  is  the  complementary  error  function.  As  an  example  of  the  Gaussian  as¬ 
sumption,  a  histogram  of  the  residual  interference  resulting  from  an  FM  interferer  is 
illustrated  in  Fig.  8  along  with  a  normal  curve.  As  can  be  seen,  the  fit  is  quite  good  for 
the  most  part.  In  fact,  for  this  example,  the  number  of  degrees  of  freedom  was  100,  and 
=  114  for  a  10%  rejection  probability.  However,  as  pointed  out  in  Section  3.2.1  in 
relation  to  the  stable  tone  for  high  signal-to-noise  ratios  at  the  output  of  the  despreader, 
performance  will  be  governed  by  the  behaviour  of  the  density  function  in  its  tail-region. 
This  behaviour  in  turn  will  also  be  a  function  of  the  operating  conditions  of  the  excisor. 


Residual  Interference 


Figure  8:  Histogram  of  residual  interference  for  the  case  when  the  interference  is  FM, 
generated  by  the  state-space  model  in  Fig.  3. 
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The  next  step  is  to  develop  5yQ  expression  for  in  Eq.  (41)  for  later 

substitution  into  Eq.  (42).  The  residual  it'terference  is,  using  Eqs.  (35)  and  (40), 

Ain  =  -  [c08(€n)  ♦  ^tPF.n]  COs(e„) 

+  l//([l/\/2nB,„  +  S„  C08(^n|n-l)]  *  COs(en)}  COs(dn) 

+  /{  [  COs(Cn)  ♦  hLPF,n]  sin(£n) 

+  l//([l/\/5nB.n +  SnCOs{4|n-l)]  ♦  ^LPF,!.)  sin(en)}sin(5„).  (43) 

Equation  (43)  shows  that  the  residual  interference  consists  of  two  terms,  one  related  to 
the  in-phase  interference  component  and  one  related  to  the  quadrature  component. 

Consider  the  in-phase  term  first.  If  the  term  [cos(c„)  *  hipp.n]  co8(c„)  is  approx¬ 
imated  by  co8*(€„)  (the  worst  case  situation  in  which  the  low  pass  filter  has  no  effect  on 
the  phase-noise  term),  and  if  c„  C  1  (see  Fig.  9  for  an  example),  and  cos(e„)  fa  1,  then 


Figure  9:  Histogram  of  phase  error  c„. 


the  in-phase  tenn  can  be  approximated  by 

Ax/,„  «  {lel  +  [1/ N/2nfl,n  +  Sn  C08(5nln-X  )]  •  ^LPF.n}  C08(fln).  (44) 

Equation  (44)  contains  two  noise  terms,  one  due  to  the  phase-noise  (a  second  order  effect), 
and  one  due  to  additive  noise  and  signal  which  have  been  low  pass  filtered. 

Now  consider  the  quadrature  term  in  Ek}.  (43).  Using  the  approximation  sin(c„)  ss 

{/€n+  ([l/v^nB,«  +  a«co8(0„i„_,)]  ♦  AtPF,n)«n}  sin(0„).  (45) 

Equation  (45)  shows  that  the  residual  interference  in  the  quadrature  direction  consists  of 
a  first  order  phase-noise  term  and  a  combination  of  the  additive  noise  and  phase- noise  (a 
second  order  effect).  Considering  Elqs.  (44)  and  (45)  one  might  intuitively  conclude  that 
the  second  order  terms  can  be  neglected  for  low  phase-noise  cases,  and  that  the  residual 
interference  is  therefore  dominated  by  the  noise  and  signal  terms  in  Eq.  (44)  and  the 
phase-noise  term  in  Ekj.  (45).  Furthermore,  for  a  narrow  bandwidth  low  pass  filter  of 
bandwidth  Bipp^  it  is  conceivable  that  the  quadrature  term  in  Eq.  (45)  could  dominate. 
This  was  in  fact  pointed  out  in  Fig.  9  of  [8]  for  the  swept  tone  case  when  little  difference 
in  performance  was  achieved  between  Bipp  =  0.05  Hz  and  0.025  Hz. 

If  it  is  assumed  that  is  weakly  correlated  with 

[l/V^nB,„  -f-  S„COs(^nln-x)]  ♦ 

in  Eqs.  (44)  and  (45),  the  residual  interference  power  in  the  in-phase  direction  is 

E{Mln]  «  {lV2)E{el}  +  P^/2  (46) 

where  is  the  power  of  the  signal-  and  noise-related  terms  in  Eq.  (44),  i.e.,  the  power 
in 

[1/V2nfl,n  +  s„  COs(5n|n-l)]  *  ^LPF.n- 

Similarly,  the  residual  interference  power  in  the  quadrature  component  is 

«  {P/2)E{el}  +  {l/2)EUl}P^ 

«  {lV2)E{el}  (47) 

where  it  has  been  assumed  that  P  Finally,  the  residual  interference  power  is, 


17 


using  Eqs.  (46)  and  (47), 

«  (IV2)E{<.1}  +  P„I2  +  (/V2)£{<) 

«  {iy2)E{cl}  +  />„/2  (48) 

where  the  term  {I^/2)E{e*]  can  be  neglected  under  low  phase-noise  conditions.  It  should 
be  noted  that  can  also  be  neglected  when  the  low  pass  filter  bandwidth  Bipf  is  small. 

Expressions  for  E{tl],  and  in  Eq.  (48)  will  now  be  determined.  If 

the  Kalman  filter  has  reached  steady  state  and  its  linear  form  is  assumed  [6],  then  the 
filtered  phase-noise  of  Eq.  (28)  is 


tn  «  <i>u,n  *  hjKal,n  (49) 

where  hKai,n  is  the  steady  state  impulse  response  of  the  optimum  filter.  Expanding 

Eq.  (49), 

=  1// [(ni,„/v^)  •  +  [s„  3in(dn)]  ♦  hKal,n] 

=  Nn/I  +  SJI  (50) 


where 


N.  =  {n',jV2)  .  Ak.u 

(51) 

Sn  ~  ('®Ti  Sm(5n)|  *  ^Kal,n- 

(52) 

Since  s„  and  from  Eq.  (25)  are  uncorrelated,  Sn  and  Nn  are  uncorrelated.  It  will 
be  assumed,  furthermore,  that  and  Sn  are  Gaussian  random  variables.  Figure  10 
illustrates  the  histogram  of  the  term  n2^J{y/2I)  in  Eq.  (51)  before  Kalman  filtering  and 
Fig.  11  illustrates  the  histogram  of  the  same  term  2d'ter  filtering.  Observe  that  in  both 
cases  the  statistics  indicate  a  definite  Gaussian  feature.  Figure  12  illustrates  the  histogram 
of  s„  sin(^n)  in  Eq.  (52)  before  filtering,  and  Fig.  13  shows  the  same  term  after  filtering. 
The  input  term  in  Fig.  13  is  far  from  Gaussian,  whereas  the  output  term  has  a  reasonable 
Gaussian  feature. 

If  the  mean-squared  values  of  the  output  phase-noise  terms  Nn/I  and  S„/I 
from  Eq.  (50)  are  combined,  they  can  be  represented  by  V'xn.*,  the  steady  state  phase- 
noise  variance  at  the  output  of  the  Kalman  filter,  the  performance  curves  of  which  were 
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Bin  Value 


Figure  10:  Histogram  of  from  Eq.  (51)  before  Kalman  filtering. 
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Filtered  Phase  N<^  (radians) 


Figure  11;  Histogram  of  ti2  ^/(V2/)  *  ^Kai,n  from  Eq.  (51). 
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Input  Phase  Noise  (radians) 


Figure  12:  Histogram  of  the  signal  term  3„  cos((?„)  from  Eq.  (52)  before  Kalman  filtering. 
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illustrated  in  Fig.  5.  With  this  Gaussian  assumption  then, 


(53) 


and 

(54) 

Finally,  to  complete  the  derivation  of  an  expression  for  in  Eq.  (41), 

an  expression  for  P^i  &om  Eq.  (48)  must  be  determined.  Recall  that  Pm  is  the  ]K>wer 
in  the  term  [l/\/2nB,„  +  s*  cos(5„|„_i)]  ♦  hcPF,n  in  Eq.  (43).  If  the  terms  n'^  „/%/2  and 
s„  cos(^n|n-i)  hrom  Eq.  (39)  are  assumed  to  be  fiat  over  the  lowpass  filter  bandwidth  with 
power  spectral  densities  Nq/2  and  0.5  respectively  (refer  to  Fig.  14  for  an  example),  then 


Figme  14:  Power  spectral  density  of  the  signal  term  s„  cos(d„|„_i)  and  noise  term  nB,„/\/2 
for  the  case  E},/Nq  =  12  dB,  (Pj  =  L  =  20). 
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Pm  ^  Blpf{NoI2  +  0.5). 


(55) 


With  these  approximations, 

E{AP}  «  )  +  Blpf{No/2  +  0.5)/2  (56) 

and 

£ 

“  Ar„/2  +  (PI2){3V,\,^  +  V,„,.)  +  +  0.5)/2 ' 

Examples  of  the  accuracy  of  the  bit  error  rate  formula  as  represented  by  Eq.  (42) 
using  Eq.  (57),  are  presented  next.  The  interference  considered  is  of  constant  envelope  and 
of  the  FM  type  generated  from  the  interference  model  in  Fig.  3.  The  interference-to-signal 
ratio  was  set  at  20  dB  and  the  bandwidth  expansion  factor  ^  was  400.  The  processing 
gain  was  £=20  (also  equal  to  Eb).  The  values  for  Eb/No  ranged  from  0  to  12  dB.  The 
modified  Kalman  algorithm  was  used  in  the  simulations.  In  addition,  the  performzmce 
curves  illustrated  in  Fig.  6  for  the  filter  form  of  the  modified  Kalman  algorithm  were  used 
to  obtiun  the  values  for  for  use  in  Eq.  (57).  Finally,  the  3  dB  bandwidth,  Bipp,  of 
a  4-pole  low  pass  Butterworth  filter  was  used  in  Eq.  (57). 

The  power  spectrum  of  the  signal,  noise,  plus  interference  is  illustrated  in  Fig.  15. 
This  was  determined  using  the  Welch  algorithm  for  two  cases,  i.e.,  EbiNo  =  0  dB  and  12 
dB. 

The  bit  error  rate  performance  is  illustrated  in  Figs.  16  to  19.  The  low  pass  filter 
bandwidth  Blpf  ranged  from  0.1  Hz  to  0,4  Hz  in  steps  of  0.1  Hz.  For  Blpf  =  0.1  and 
0.2  Hz  (Figs.  16  and  17),  there  is  reasonable  agreement  between  theory  and  the  computer 
simulations.  The  difference  starts  to  manifest  itself  for  larger  values  of  Blpf^  i.e.,  0.3  and 
0.4  Hz  (Figs.  18  and  19).  This  is  to  be  expected,  as  it  was  pointed  out  in  the  development 
of  E)q.  (57)  that  the  correlation  between  the  low  pass  filtered  signal  and  noise  terms,  (i.e., 
the  term 

{[nB,J y/2  +  SnC08(d„|„_i)]  ♦  hi,PF,n}  C08(d„|„_,) 

in  Eqs.  (44)  and  (45))  and  s„+n„  would  increase  with  increasing  B^pp.  This  phenomenon 
is  illustrated  in  Figs.  20  and  21  for  Blpf  =  0.1  and  0.4  Hz,  respectively. 
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Frequency  (Hz) 


Figure  15:  Power  spectral  density  of  s„  +  ««  +  *n  for  E^INq  =  0  dB  and  ^'j/A'o  =  12  dB 


EhINo  (dB) 

Figure  16:  Comparison  of  the  bit  error  rate  performance  of  the  modified  Kalman  filter 
algorithm  with  Eq.  (42)  for  the  case  of  =  0.1  Hz. 
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Figxire  17:  Comparison  of  the  bit  error  rate  performance  of  the  modified  Kalman  filter 
algorithm  with  Eq.  (42)  for  the  case  of  Blpf  =  0-2  Hz. 


E,INo  (dB) 


Figure  18:  Comparison  of  the  bit  error  rate  performance  of  the  modified  Kalman  filter 
algorithm  with  Eq.  (42)  for  the  case  of  B^pp  —  0.3  Hz. 


£j/jVo  (dB) 

Figure  19:  Comparison  of  the  bit  error  rate  performance  of  the  modified  Kalman  filter 
algorithm  with  Eq.  (42)  for  the  case  of  B[,pp  =  0.4  Hz. 
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Figure  20  shows  very  little  distortion  of  the  signal  at  the  output  of  the  excisor, 
i.e., 

s„  -  {[Sn  CO8(0„|„_i)]  ♦  hi,PF,n}  COs(^„|„_i  ) 
for  Blpf  =  0.1  Hz  before  the  PN  correlator  in  Fig.  1.  The  curve  referring  to 

{[SnC08((?n|n_i)]  ♦  hLPF,n}  C08(5„|„_i), 


is  approximately  10  dB  lower  in  power  in  the  vicinity  of  0  Hz.  Figure  21 ,  which  corresponds 
to  Blpf  =  0.4  Hz,  shows  some  distortion  of  the  signal.  In  fact,  the  difference  in  power 
was  20%  or  1  dB.  Similar  distortion  phenomena  also  occurred  with  the  noise  term 

{«B.n  ♦  ^LPF,n]  C08(d„|n_i ) 

for  the  two  lowpass  filter  bandwidths,  although  these  are  not  illustrated  here. 

Finally,  Fig.  22  shows  an  example  of  the  power  spectral  density  of  the  interfer¬ 
ence,  t'n,  interference  estimate, 

A 

{{/C0s(£n)]  ♦  htPF.n}  COs(0„|„.l) 

(which  excludes  the  signal  and  noise  term  effects),  and  the  residual  interference, 

A 

~  {[/cOs(Cn)]  ♦  hLPF,n}  CO8(0„|„_i), 

for  Ei,/Nq  =  12  dB  and  Blpf  =  0.1  Hz.  The  degree  of  interference  suppression  for  this 
example  was  17.38  dB,  and  the  mean-squared  value  of  the  residual  inteifcieuce  was  1.83. 
For  Blpf  =  0.4  Hz,  the  degree  of  interference  suppression  and  mean-squared  value  of 
the  residual  interference  were  almost  identical  to  that  for  Blpf  =  0.1  Hz,  i.e.,  17.39  dB 
and  1.82.  This  leads  one  to  the  conclusion  that  residual  interference  was  not  a  factor  in 
the  reasons  for  the  difference  between  the  bit  error  rate  curves  of  Blpf  =  0.4  Hz,  but 
that  it  was  due  to  other  factors  such  as  signal  and  noise  distortion  and  the  violation  of 
several  assumptions  which  were  made  in  the  development  of  the  bit  error  rate  formula  in 
Eq.  (42),  the  primary  one  being  that  Blpf  be  significantly  less  than  0.50  Hz. 
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A 

Figure  20:  Power  spectral  density  of  signal  Sn,  the  filtered  signal  term  {[5„  cos(5„l„_i)]  ♦ 
hLPF,n}  co8(0„|„_i)  and  the  residual  signal  term  -  {[sn  cos(^„|n-i)]  *  ^LPF.n}  cos(^„i„_,). 
The  bandwidth  of  the  lowpass  filter  is  Bcpf  —  0.1  Hz. 
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10» 


Figure  21:  Power  spectral  density  of  signal  Sn>  the  filtered  signal  term  {fs„  cos(^„(„_i)]  * 

A  A  * 

h£,PF,n}cos(dn|n-i),  and  the rcsidual signal  tenns„-{[snCos(0n|n-i)]*hi,PF,n}cos{0n|n-i). 
The  bandwidth  of  the  lowpass  filter  is  Blpf  —  0.4  Hz. 
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Figure  22:  Power  spectral  density  of  interference  t„,  the  interference  estimate  {[/  cos(cn)]  * 
co8(0n|n-i),  and  the  rcsidual  termin-{[/co8(c„)|*/»LPr,n}cos(5„l„_i).  The  band¬ 
width  of  the  lowpass  filter  is  Blpf  =  0.1  Hz. 


6.0  CONCLUSIONS 

A  bit  error  rate  expression  for  the  Kalman  filter  excisor  has  been  developed  based  on 
several  assumptions  which  were  elaborated  upon.  The  interference  in  this  case  was  of  the 
FM-type  generated  by  the  state  space  model  on  which  the  Kalman  filter  algorithm  was 
derived.  It  was  shown  and  illustrated  that  theory  agrees  well  with  simulation  as  long  as  the 
bandwidth  Bcpf  of  the  low  p^s  filter  used  to  determine  the  envelope  of  the  interference 
is  much  less  than  the  sampling  rate.  It  was  edso  shown  that  the  interference  estimate 
at  the  input  to  the  interference  canceller  will  have  significant  cross-correlation  with  the 
signal,  leading  to  some  cancellation  of  the  spread  spectrum  signal,  if  the  bandwidth  jBlpf 
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is  large  relative  to  the  sampling  rate. 

For  constant  envelope  interference,  or  interferers  whose  envelopes  vaury  slowly 
with  time,  the  approach  taken  in  this  report  will  work  well.  However,  for  narrowband 
Gaussian  noise  it  will  not,  this  being  due  to  the  sudden  phaise  changes  when  the  envelope 
goes  to  zero  and  the  resulting  temporary  loss  of  lock  of  the  Kalman  filter.  The  effect 
would  be  an  increase  in  residual  interference  which  could  not  be  modelled  accurately  as 
Gaussian.  More  work  in  this  area  is  required. 
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